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Abstract 



We evaluate the continuum limit of the valence (quenched) approxima- 
tion to the mass of the lightest scalar quarkonium state for a range of dif- 



^^ I ferent quark masses and to the mixing energy between these states and the 

_i. ' lightest scalar glueball. Our results support the interpretation of /o(1710) as 

composed mainly of the lightest scalar glueball. 
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I. INTRODUCTION 

Evidence that /o(1710) is composed mainly of the hghtest scalar glueball is now given 
by two different sets of numerical determinations of QCD predictions using the theory's 
lattice formulation in the valence (quenched) approximation. A calculation on GFll [^ of 
the width for the lightest scalar glueball to decay to all possible pseudoscalar pairs, on a 
16^ X 24 lattice with /3 of 5.7, corresponding to a lattice spacing a of 0.140(4) fm, gives 
108(29) MeV. This number combined with any reasonable guesses for the effect of finite 
lattice spacing, finite lattice volume, and the remaining width to multibody states yields a 
total width small enough for the lightest scalar glueball to be seen easily in experiment. For 
the infinite volume continuum limit of the lightest scalar glueball mass, a reanalysis of 
a calculation on GFll 0, using 25000 to 30000 gauge configurations, gives 1648(58) MeV. 
An independent calculation by the UKQCD-Wuppertal [Q collaboration, using 1000 to 3000 
gauge configurations, when extrapolated to the continuum limit according to Refs. |Q,|^ 
yields 1567(88) MeV. A more recent calculation using an improved action ^ gives 1730(94) 
MeV. The three results combined become 1656(47) MeV. A phenomenological model of the 
glueball spectrum which supports this prediction is discussed in Ref. |0]. 

Among established resonances with the quantum numbers to be a scalar glueball, all are 
clearly inconsistent with the mass calculations except /o(1710) and /o(1500). Between these 
two, /o(1710) is favored by the mass result with largest statistics, by the combined result, 
and by the expectation [^ that the valence approximation will lead to an underestimate of 
the scalar glueball's mass. Refs. |jl|J^ interpret /o(1500) as dominantly composed of strange- 
antistrange, ss, scalar quarkonium. A possible objection to this interpretation, however, is 
that /o(1500) apparently does not decay mainly to states containing an s and an s quark 0. 
In part for this reason, Ref. 0] interprets /o(1500) as composed mainly of the lightest scalar 
glueball and /o(1710) as largely ss scalar quarkonium. A second objection is that while the 
Hamiltonian of full QCD couples quarkonium and glueballs, so that physical states should be 
linear combinations of both, mixing is not treated quantitatively in Ref. |]I[ . In the extreme, 
mixing could lead to /o(1710) and /o(1500) each half glueball and half quarkonium. 

Using the valence approximation for a fixed lattice period L of about 1.6 fm and a range of 
different values of quark mass, we have now calculated the continuum limit of the mass of the 
lightest scalar qq states and the continuum limit of the mixing energy between these states 
and the lightest scalar glueball. Our calculations have been done with four different choices 
of lattice spacing. Continuum predictions are found by extrapolation of results obtained 
from the three smallest values of lattice spacing. For the two choices of lattice spacing we 
have also calculated scalar qq masses on lattices with L of about 2.3 fm, and for one choice of 
lattice spacing we have found scalar quarkonium-glueball mixing energies on a lattice with 
L of about 2.3 fm. Preliminary versions of this work are reported in Refs. P,0,Q. 

Our results provide answers to the objections to the interpretation of /o(1710) as largely 
the lightest scalar glueball. For the valence approximation to the infinite volume continuum 
limit of the ss scalar mass we find a value significantly below the valence approximation 
scalar glueball mass. This prediction makes improbable, in our opinion, the identifica- 
tion |l^ of /o(1500) as primarily a glueball and /o(1710) as primarily s's quarkonium. Our 
calculation of glueball-quarkonium mixing energy, combined with the simplification of con- 
sidering mixing only among the lightest discrete isosinglet scalar states, then yields a mixed 



/o(1710) which is 73.8(9.5)% glueball and a mixed /o(1500) which is 98.4(1.4)% quarkonium, 
mainly ss. The glueball amplitude which leaks from /o(1710) goes almost entirely to the 
state /o(1390), which remains mainly nn, normal-antinormal, the abbreviation we adopt for 
{uu + dd)/\/2. We find also that /o(1500) acquires an nn amplitude with sign opposite to 
its ss component suppressing, by interference, the state's decay to KK final states. Assum- 
ing SU(3) flavor symmetry before mixing for the decay couplings of scalar quarkonium to 
pairs of pseudoscalars, the KK decay rate of /o(1500) is suppressed by a factor of 0.39(16) 
in comparison to the rate of an unmixed ss scalar. This suppression is consistent, within 
uncertainties, with the experimentally observed suppression. 

It perhaps is useful to discuss briefly at this point a proposed calculation of mixing be- 
tween valence approximation quarkonium and glueball states through common decay chan- 
nels [0 which forms the basis for additional objections to the identification of /o(1710) as 
primarily a glueball. A detailed examination of problems with the calculation of Ref. |T^ 
appears in Ref. [13 ]. One defect of the work in Ref. |13] is the omission of quarkonium 



to glueball transitions by direct annihilation of the quarkonium's quark and antiquark into 
chromoelectric field. Direct annihilation is the leading valence approximation contribution 
to mixing and is evaluated in the present paper. On the other hand, the transitions through 
two-pseudoscalar intermediate states which Ref. |]13[ considers include an extra closed quark 
loop in addition to the quark paths of the direct quark-antiquark annihilation process. Thus 
according to a systematic scheme for evaluating all quark loop corrections to the valence 
approximation |1^ , the decay channel mixing calculation is part of the one-quark-loop cor- 



rection to the direct annihilation mixing amplitude. As shown in detail in Ref. [Q, however, 
the corrections to the direct annihilation amplitude must include also a counterterm pro- 
portional to the pure gauge action. This counterterm is required to compensate for the 
shift between the screened effective gauge coupling used in the valence approximation and 
the unscreened bare coupling of full QCD. The counterterm is entirely absent from the 
calculation of Ref. [O]. As a consequence of this omission and of the omission of the di- 



rect quark-antiquark annihilation term, we believe the mixing calculation of Ref. ||T3[ is not 
correct. 

The remainder of this paper is organized as follows. In Section || we describe the Monte 
Carlo ensembles of gauge field configurations we use. In Section |TTl| we present the calculation 
of scalar quarkonium masses. In Section ^V\ we describe a glueball mass calculation. In 
Section we present a calculation of quarkonium-glueball mixing energy. In Section ^ 
we consider the physical mixed glueball and quarkonium states. Finally, Section |V11| briefly 
examines consequences of quarkonium-glueball mixing for glueball decay. 

II. MONTE CARLO ENSEMBLES 

Our calculations, using Wilson fermions and the plaquette action, were done for four 
choices of (3 with two different lattice sizes at each of the two smallest (3 giving a total of six 
combinations of (3 and lattice structure. These are listed in Table |I[ For each combination 
of P and lattice structure, calculations were typically done with five different choices of k. 
These are listed in Table |I| along with the corresponding Monte Carlo ensemble sizes. In all 
cases, a sufficient number of updating sweeps was skipped between successive Monte Carlo 
configuration to leave no statistically significant correlations between successive pairs. The 



ensemble of 599 configurations used for two hopping constant choices at /3 of 6.4 is a subset 
of the 1003 configuration ensemble used for the three other k at /? of 6.4. At (3 of 5.70 on a 
lattice 16^ x 24, the 3870 member ensemble at k, of 0.1600 is a subset of the 12186 member 
ensemble at k, of 0.1650, which has no overlap with the 1972 member ensemble used at k, of 
0.1625 and 0.1650. For all other entries in Table |I|, all n values share a single ensemble of 
gauge field configurations. 

From the smallest to largest (3, the lattice spacing varies by nearly a factor of 2.7. The 
smaller lattices with (3 of 5.70 and 5.93, and the lattices with /3 of 6.17 and 6.40 have 
nearly the same periods in the two (or three) equal space directions and thereby permit 
extrapolations to zero lattice spacing with nearly constant physical volume. 

The values of lattice spacing and lattice period in Table | and conversions from lattice to 
physical units in the remainder of this article, are determined MM from the exact solution 



to the two-loop zero-flavor Callan-Symanzik equation for A.jj-a with A^^ of 234.9(6.2) 



MS MS 



MeV determined from the continuum limit of (Ajj-a) / {nipa) in Ref. |jT5[. For P from 5.70 



to 6.17, the ratio (Ai^a)/(mpa) in Ref. [|1^ was found to be constant within statistical 
errors, thus our results are, within errors, almost certainly the same as those we would have 
obtained by converting to physical units using values of rripa. We chose to convert using 
Aj^, however, since Ref. [|l^ did not find TUpU at (3 of 6.40, which would be needed for our 
present calculations. 

III. QUARKONIUM MASSES 

For each ensemble of gauge fields, with two exceptions at (3 of 5.70, we evaluated cor- 
relation functions using smeared Coulomb gauge quark and antiquark fields incorporating 



random sources following Ref. |T^ . From these fields we constructed pseudoscalar and scalar 
quarkonium propagators. Averaged over the random sources, the propagators we calculated 
become 

Cjfit) = Y.<f,ix,t)Mo,o)>, (1) 

X 

where / is either tt or a and 7r„(x, t) and o"j,(x, t) are, respectively, the smeared pseudoscalar 
and scalar operators of Ref. |T^ with smearing size v. At (3 of 5.70 on a lattice 16^ x 24, for 
the 3870 member gauge ensemble with k of 0.1600 and for the 12186 member ensemble at k 
of 0.1650, the propagators of Eq. (|l]) were evaluated directly without use of random sources. 
Evidence in Ref. [O suggests that for equal statistical uncertainties, propagators found 



using random sources require about half the computer time needed for a direct calculation. 
The values used for r for each (3 and lattice are listed in Table |T|. 

For sufficiently large values of t and the lattice time period T, C.j,T,(t) and Ccraif) are 
expected to approach the asymptotic form 

Cffit) — > Zf[exp{—mfat) + expimjat — mfoT)], (2) 

where / can be either vr or a. Fitting the large t behavior of C,r7r(^) and C„„{t) to Eq. (|^) we 
obtained the masses, in lattice units, m.,^a and maa, and the field strength renormalization 
constants Z^, and Z„. 



Fitting CT,T,{t) and Caaif) to Eq. (Q) at pairs of neighboring time slices t and t + 1 gives 
the effective masses 171^^(1) and m„{t), which at large t approach m^ and m^, respectively. 
To determine ?7i^, Z^r, mg-, and Z^-, we began by examining effective mass graphs for a range 
of g, in Eq. ([^), to find smearing sizes for which m.,t{t) and m„{t) show clear evidence of 
approaching constants at large t. In all but one case we found satisfactory effective mass 
plateaus with q the same as r of Table |T|. For C^^ at j3 of 6.17 a smearing size of 9.0 was 
used for q. Typical effective mass graphs are shown in Figures |I] - |T^. 

Trial time intervals on which to fit C.,rT^{t) and Cfj„{t) to Eq. (Q) were chosen from effective 
mass graphs by eliminating large values of t with large statistical uncertainties in effective 
masses and eliminating small t at which effective masses have clearly not yet reached the 
large t plateau. Fits were then made to Eq. (Q) on all subintervals of 3 or more consecutive 
t within the trial range. The fit for each interval was chosen to minimize the x^ taking 
into account all correlations among the fitted data. Correlations were determined by the 
bootstrap method. The final fitting interval for each propagator was chosen to be the interval 
with the smallest x^ P^r degree of freedom. 

Final fitting intervals and fitted masses are shown by solid lines in Figures ^ - 11^. Dashed 
lines extend the solid lines toward smaller times to display the approach of effective masses 
to the final fitted masses. 

Tables |T^ - pCVll list the final pseudoscalar and scalar masses obtained. The statistical 



uncertainties for the masses in these tables, and in all other Monte Carlo results in this 
article, are determined by the bootstrap method. 



For (3 of 5.70, 5.93, 6.17 and 6.4, Figures 0, |T8|, |19|, and ^ respectively, show the 
pseudoscalar mass squared ml as a function of 1/k. The solid line in each figure shows a fit 
of ml to a quadratic function of 1/k, used to determine the strange quark hopping constant 
Kg at which 

ml = 2Ml - Ml (3) 

where Mk and M,r are the observed neutral kaon and pion masses, respectively. The 
quadratic fits in I/k were used also to determine the critical hopping constant Kent at 
which m^r is zero. Although the determination of Kf-nt depends on extrapolation of each fit 
beyond the kappa interval in which we have data, the determination of Kg does not and uses 
the fits only to interpolate between measurements. From Kcru we define the quark mass for 
each K to be 



Values of Kg and Kctu are given in Table pCVlll. 



For the two lattices with (3 of 5.93, Figure ^ shows the scalar quarkonium mass as 
a function of quark mass ^a. The solid lines in Figure ^ are fits of the scalar mass to 
quadratic functions of quark mass. The scalar masses found by interpolation to the strange 
quark mass are also indicated. As shown by the figure, for the lattice 16^ x 14 x 20 with 
L of 1.54(4) fm the scalar mass as a function of quark mass flattens out as quark mass is 
lowered toward the strange quark mass and then appears to begin to rise as the quark mass 
is decreased still further. This feature is absent from the data at (3 of 5.93 for the lattice 
24^ with L of 2.31(6) fm and is thus a finite- volume artifact. It is present in the data at (3 



of 5.70 with L of 1.68(5) fm, at p of 6.17 with L of 1.74(5) fm, and at p of 6.40 with L of 



1.66(5) fm shown in Figures ^, 23 and |2J, respectively. It is absent, however, in the data 



at /3 of 5.70 with L of 2.24(7) fm shown in Figure |2^. Values of the scalar quarkonium mass 



interpolated to the strange quark mass are given in Table PCIX 

The pseudoscalar mass squared m^ shown in Figures |T^ - ^ is nearly a linear function 
of I/k for all (3 and lattice periods. The difference in rriTra between the two lattice at (3 of 
5.70 and between the two lattice at (3 of 5.93 is in all cases less than 0.5%. The anomaly 
in quark mass dependence of the scalar mass for L of 1.6 fm, shown in Figure ^, is absent 
from the quark mass dependence of the pseudoscalar mass for this value of L. 

For L near 1.6 fm. Figure ^ shows the ss scalar mass in units of Ay^ as a function of 

lattice spacing in units of 1/A^n4r. A linear extrapolation of the mass to zero lattice spacing 
gives 1322(42) MeV, far below our valence approximation infinite volume continuum glueball 
mass of 1648(58) MeV. For the ratio of the ss mass to the infinite volume continuum limit 
of the scalar glueball mass we obtain 0.802(24). Figure ^ shows also values of the ss scalar 
mass at (3 of 5.70 and 5.93 with L of 2.24(7) and 2.31(6) fm, respectively. The ss mass 
with L near 2.3 fm lies below the 1.6 fm result for both values of lattice spacing. Thus the 
infinite volume continuum ss mass should lie below 1322(42) MeV. We believe our data make 
improbable the interpretation of /o(1500) as mainly composed of the lightest scalar glueball 
with /o(1710) consisting mainly of ss scalar quarkonium. For comparison with our data. 
Figure ^ shows the valence approximation value for the infinite volume continuum limit of 
the scalar glueball mass and the observed value of the mass of /o(1500) and of the mass of 
/o(1710) The uncertainties shown in the observed masses in units of A^^ arise mainly from 
the uncertainty in A^;^. 

IV. GLUEBALL MASS 

In preparation for a calculation of quarkonium-glueball mixing energy, from each gauge 
ensemble we also constructed scalar glueball operators. On the gauge ensembles at /5 of 5.70, 
we evaluated smeared Coulomb gauge scalar glueball operators and at all larger /3 smeared 
gauge invariant scalar glueball operators. The operators we used are discussed in Ref. 0. 
The correlation function constructed from these is 

Cgg{t) = -Y,[< 9{x,t)g{y,0) >- < g{x,t) >< g{y,0) >], 

xy 

(5) 

where g{x,t) is the smeared scalar glueball operator and V is the space direction lattice 
volume. 

Fitting the the large t behavior of Cgg(t) to Eq. (0) for / chosen to be f?, we obtain the 
glueball mass rriga and field strength renormalization constant Zg. A detailed discussion 
of calculations of rUga and Zg for the same (3 and nearly the same lattice sizes considered 
here, but with much larger Monte Carlo ensemble sizes, is presented in Ref. [§]. Using the 
calculation of Ref. pi to guide the choice of smearing parameters and time intervals to be 



fit, we applied the fitting procedure of Section |ITT[ Smearing parameters we found to be 



satisfactory are given in Table pT[ Effective mass graphs are shown in Figures |26| - |30. 



Fitted masses, fitted time intervals, y^ per degree of freedom of each fit, and corresponding 
lattice sizes and fitted masses from Ref. are given in Table PCX| . For the Monte Carlo 
ensemble with /5 of 5.70 on a lattice 16^ x 24 the fits of Cgg{t) to Eq. (Q) yielded either large 
X^ or large statistical errors, therefore no results are given in Table PCX| for this case. 

V. MIXING ENERGY 

To determine scalar quarkonium-glueball mixing energies, we evaluated the correlation 



between the scalar quarkonium operators of Section |TT| and the glueball operators of Sec- 
tion 1^. For scalar quarkonium operators not containing random variables and for the ran- 
dom operators when averaged over random variables, the correlation function we calculated 
becomes 

CgAt) = E[< 9{x,t)a{0,0) >-< g{x,t) >< a(0,0) > . (6) 

X 

The smearing parameters for quark and glueball fields, as before, are listed in Table |T^. For 
large t and time period T, the asymptotic behavior of Cg„{t) for m„ close to nig is 



Cga{t) -^ J ZgZsEa'^[exp{—mga\t — t'\) + exp{mga\t — t'\ — rrigaT)] x 

[exp{—ma-a\t'\) + exp{m„a\t'\ — m^aT)]. (7) 



Fitting Cg„{t) to Eq. (0) using m„, Z^, nig and Zg from Sections [IT^ and |^, we found the 
glueball-quarkonium mixing energy in lattice units Ea. To choose the t range over which to 
fit Cgait) to Eq. d^), it is convenient to define an effective mixing energy E{t) by fitting Cga(t) 
to Eq. (0) solely at t. Typical data for E{t) is shown in Figures ^- ^ Trial time intervals 
on which to fit Cga{t) to Eq. (^ were chosen from graphs of E{t), following the fitting 
procedure of Section |T|, by eliminating large values of t with large statistical uncertainties 
in E{t) and eliminating small t at which E{t) has clearly not yet reached a large t plateau. 
Fits with minimal correlated x^ were then made to Eq. (^ on all subintervals of 2 or more 
consecutive t within the trial range. The final fitting interval for each propagator was chosen 
to give the smallest x^ per degree of freedom. 



Final mixing energy values are given in Tables pCXl| - pCXV|. A few of the combinations of 



/?, K and lattice size appearing in Table |T| are missing from Tables |XX1] - PCXV| . No results 
are given for /3 of 5.7 on the lattice 16^ x 24 since, as mentioned in Section fV|, we were 
unable to obtain stable values for nig and Zg for this data set. We also give no results for j3 
of 5.7 and n of 0.1650 on the lattice 12^ x 10 x 24, for which the scalar quarkonium fit was 
poor, and no results for (3 of 6.4 and n of 0.1485 and 0.1488 on 32 x 28 x 40, for which the 
Monte Carlo ensembles were too small to give reliable values of Cga{t). 



Figure |3^ shows the quarkonium-glueball mixing energy as a function of quark mass for 
the two different lattices with (3 of 5.93. For neither lattice does there appear to be any 
sign of the anomalous quark mass dependence found in Figure ^ The mixing energies at 
different quark masses turn out to be highly correlated and depend quite linearly on quark 
mass. Figures ^, ^ and ^ show mixing energy as a function of quark mass for j3 of 5.70, 
6.17 and 6.40, respectively. For these values of j3 the mixing energy also shows no sign of the 



anomalous quark mass dependence exhibited by the scalar quarkonium mass. The nearly 
linear dependence of Figure ^ is repeated. Thus it appears that the mixing energy can be 
extrapolated reliably down to the normal quark mass //„, defined to be the quark mass at 
which 771^ becomes M^. 

Table pCXVl] gives values of the mixing energy interpolated to the strange quark mass 
fis, extrapolated down to the normal quark mass fin, and of the ratio of these two energies. 
For the data at /3 of 5.93, the ratio changes by less than 3% from L of 1.54(4) fm to L of 
2.31(6) fm, a difference consistent with the statistical error. Thus the ratio has at most small 
volume dependence and seems already to be near its infinite volume limit with L around 1.6 
fm. 



Figure |4^ shows linear extrapolations to zero lattice spacing of quarkonium-glueball 
mixing energy at the strange quark mass -E(/is) and of the ratio E (fin) /E (fig)- The zero 
lattice spacing prediction E{fis) is 43(31) MeV and of £"(//„ )/i?(/is) is 1.198(72). 

VI. MIXED PHYSICAL STATES 

We now combine our infinite volume continuum value for E{^n)/E{^s) with a simplified 
treatment of the mixing among valence approximation glueball and quarkonium states which 
arises in full QCD from quark-antiquark annihilation. The simplification we introduce is to 
permit mixing only between the lightest scalar glueball and the lowest lying discrete quarko- 
nium states. We ignore mixing between the lightest glueball and excited quarkonium states 
or multiquark continuum states, and we ignore mixing between the lightest quarkonium 
states and excited glueball states or continuum states containing both quarks and glueballs. 

Excited quarkonium and glueball states and states containing both quarks and glueballs 
are expected to be high enough in mass that their effect on the lowest lying states will 
be much smaller than the effect of mixing of the lowest lying states with each other. On 
the other hand, as mentioned earlier, according to the systematic version of the valence 



approximation described in Ref. [|1J], the additional feedback into mixing among the lowest 
discrete quarkonium and glueball states arising as a consequence of the coupling, omitted 
from our simplified mixing, of the lowest glueball and scalar quarkonium states to continuum 
multi-meson states is a quark loop correction to the direct glueball-quarkonium mixing 
amplitude which our simplified mixing includes. For low energy QCD properties there is 
a reasonable amount of phenomenological evidence that such quark loop corrections are 
relatively small. 

The structure of the Hamiltonian coupling together the scalar glueball, the scalar ss and 
the scalar nn isosinglet becomes 



m 



g 



E{fi,) V2rE{fis 
E{fis) m^ifJ's) 

V2rE{fis) m„{nn). 

Here r is the ratio E{^n)/E{fis) which we found to be 1.198(72), and nig, m„{^s) and m^{^n) 
are, respectively, the glueball mass, the ss quarkonium mass and the nn quarkonium mass 
before mixing. 

The three unmixed mass parameters we take as unknowns. We will also treat E{fis) as an 
unknown since the fractional error bar on our measured value is large. The four unknown 
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parameters can now be determined from four observed masses. To leading order in the 
valence approximation, with valence quark- ant iquark annihilation turned off, corresponding 
isotriplet and isosinglet states composed of u and d quarks will be degenerate. For the 
scalar meson multiplet, the isotriplet {uu — dd)/\/2 state has a mass reported by the Crystal 
Barrel collaboration to be 1470(25) MeV §]. Thus we take m^{fXn) to be 1470(25) MeV. In 
addition, the Crystal Barrel collaboration finds an isosinglet mass of 1390(30) MeV from 
one recent analysis and 1380(40) MeV 0] from another. Mark III finds 1430(40) MeV []T7|. 



We take the mass of the physical mixed state with largest contribution coming from nn 
to be 1404(24) MeV, the weighted average of 1390(30) MeV and 1430(40) MeV. The mass 
of the physical mixed states with the largest contributions from ss we take as the mass of 
/o(1500), for which the Particle Data Group's averaged value is 1505(9) MeV. The mass 
of the physical mixed state with the largest contributions from the glueball we take as the 
Particle Data group's averaged mass of /o(1710), 1697(4) MeV. 

Adjusting the parameters in the matrix to give the physical eigenvalues we just specified, 
rUg becomes 1622(29) MeV, m^{fis) becomes 1514(11) MeV, and E{fis) becomes 64(13) MeV, 
with error bars including the uncertainties in the four input physical masses. The unmixed 
iTig is consistent with the world average valence approximation glueball mass 1656(47) MeV, 
E{fis) is consistent, within large errors, with our measured value of 43(31) MeV, and m„{ns) 
is about 13% above the valence approximation value 1322(42) MeV for lattice period 1.6 
fm. This 13% gap is comparable to the largest disagreement, about 10%, found between 
the valence approximation and experimental values for the masses of light hadrons. As 
expected from the discussion of Ref. p|], the valence approximation value lies below the 
number obtained from experiment. 

For the three physical eigenvectors we obtain 

|/o(1710) > = 0.859(54)1^ > +0.302(52)|ss > +0.413(87) |nrl >, 

|/o(1500) > = -0.128(52)1^ > +0.908(37)|ss > -0.399(113)|nrr >, (8) 

|/o(1390) > = -0.495(118)1^ > +0.290(91)|ss > +0.819(89) |nrr > . 

The mixed /o(1710) has a glueball content of 73.8(9.5)%, the mixed /o(1500) has a glueball 
content of 1.6(1.4)% and the mixed /o(1390) has a glueball content of 24.5(10.7)%. Since, 
as well known, the partial width r(J/\E' — ^ 7 + /i) is a measure of the size of the gluon 
component in the wave function of hadron h, our results imply that r(J/\l/ — > 7 + /o(1710)) 
should be significantly larger than r(J/^ ^ 7 + /o(1390)) and T{J/^ -> 7 + /o(1390)) 
should be significantly larger than r(J/^ -^ 7 + /o(1500)). These predictions are supported 
by a recent reanalysis of Mark III data []T^. In addition, in the state vector for /o(1500), 
the relative negative sign between the ss and nn components will lead, by interference, to 
a suppression of the partial width for this state to decay to KK. Assuming SU(3) flavor 
symmetry for the two pseudoscalar decay coupling of the scalar quarkonium states, the 
total KK rate for /o(1500) is suppressed by a factor of 0.39(16) in comparison to the KK 
rate for an unmixed ss state. This suppression is consistent, within uncertainties with the 
experimentally observed suppression. 



VII. GLUEBALL DECAY COUPLING 

We now consider briefly the contribution to scalar glueball decay to pseudoscalar quarko- 
nium pairs arising from quarkonium-glueball mixing. 

In Ref. [m a calculation of scalar decay to pseudoscalar quarkonium pairs was done on 
a spatial lattice of 16^ at (3 of 5.70 and k of 0.1650 and 0.1675. For these parameters, the 
results of Section |T| imply the lightest scalar quarkonium state is significantly heavier than 
the lightest scalar glueball. It is not hard to show that in this circumstance, the valence 
approximation decay calculation includes, to first order in the quarkonium-glueball mixing 
energy, the contribution arising from mixing of the scalar glueball with scalar quarkonium. 
This first order contribution is 

771 

AA(0 ^ TT + Tt), = A(cr ^ TT + Tt), (9) 

nifj — TUg 

where, as before, a is the lightest scalar quark-antiquark state and tt is the lightest pseu- 
doscalar quark-antiquark state all with a single common value of k. 

Although we do not have values for X{a ^ tt + tt) at /3 of 5.70, a rough estimate of the 
order of magnitude of AA(cr ^ tt + tt) can be made by taking X{a —>■ tt + tt) from experiment. 
Assuming SU(3) flavor symmetry for scalar quarkonium decay couplings, the observed decay 
width of the scalar i^*(1430) yields A(cr — *> vr + vr) of about 8 GeV. Combining this number 
with Ea of about 0.2, and mrj — mg of about 0.3, we get AX{g ^ tt + tt) of about 5 GeV. The 
A((7 ^ TT + tt) found in Ref. DH range from about 1.5 to 3 GeV. It thus highly probable that 
the glueball decay couplings of Ref. |1[] include significant contributions from mixing of the 
scalar glueball with scalar quarkonium. It appears possible that the decay couplings may 
arise entirely from the mixing contribution. A lattice calculation of X{a ^ tt + tt) would 
confirm or refute this possibility. If glueball decay were found at (3 of 5.70 to occur entirely 
through mixing, a reasonable guess would be that this is also the decay mechanism in the 
real world. 
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TABLES 



/? 


lattice 


a (fm) 


period (fm) 


5.70 


122 X 10 X 24 


0.140(4) 


1.68(5) 


5.70 


16^ X 24 


0.140(4) 


2.24(6) 


5.93 


16^ X 14 X 20 


0.0961(25) 


1.54(4) 


5.93 


24^ 


0.0961(25) 


2.31(6) 


6.17 


242 X 20 X 32 


0.0694(18) 


1.74(5) 


6.40 


32^ X 28 X 40 


0.0519(14) 


1.66(5) 



TABLE I. For each (3 and lattice structure, the corresponding lattice spacing and lattice period 
in the two (or three) equal space directions in fernii. 



/? 


lattice 


K 


ensemble size 


5.70 


122 X 10 X 24 


0.1600, 0.1613, 0.1625, 0.1638 


2749 


5.70 


16^ X 24 


0.1600 
0.1625 
0.1650 


3870 

1972 

1972, 12186 


5.93 


16^ X 14 X 20 


0.1539, 0.1546, 0.1554, 0.1562, 0.1567 


2328 


5.93 


244 


0.1539, 0.1554, 0.1567 


1733 


6.17 


242 X 20 X 32 


0.1508, 0.1523, 0.1516, 0.1520, 0.1524 


1000 


6.40 


322 X 28 X 40 


0.1485, 0.1488 
0.1491, 0.1494, 0.1497 


599 
1003 




TABLE IL Monte Carlo ensemble size for each (3, lattice structure. 


and K. 


/? 


lattice 


quark smearing N 


S C 


5.70 


122 X 10 X 24 


2.0 


1 


5.70 


16^ X 24 


2.0 


1 


5.93 


162 X 14 X 20 


3.0 7 


6 


5.93 


24^ 


3.0 7 


6 


6.17 


242 X 20 X 32 


4.5 7 


7 


6.40 


322 X 28 X 40 


6.0 8 


9 



TABLE in. Quark smearing parameters, glueball gauge invariant smearing parameters N and 
S, and glueball Coulomb gauge smearing parameter C. 
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mass 



t range 



xVd.o.f. 



0.1600 
0.1613 
0.1625 
0.1638 
0.1650 



0.6884(8) 
0.6330(8) 
0.5795(8) 
0.5176(10) 
0.4549(11) 



7 



9 



10 
11 
11 
10 
11 



0.03 
0.06 
0.08 
0.01 
0.00 



TABLE IV. For /3 of 5.70 on a lattice 12^ x 10 x 24, for each value of k, fitted pseudoscalar 
meson masses, time range of fit, and fit's x^ per degree of freedom. 



mass 



t range 



x'/d.o.f. 



0.1625 
0.1650 



0.5795(4) 
0.4560(5) 



7- 10 
7- 10 



0.37 
0.25 



TABLE V. For (3 of 5.70 on a lattice 16^ x 24, for each value of k, fitted pseudoscalar meson 
masses, time range of fit, and fit's x^ per degree of freedom. 



mass 



t range 



xVd.o.f. 



0.1600 
0.1650 



0.6888(5) 
0.4572(3) 



6 

7 



9 



0.00 
0.08 



TABLE VI. For (3 of 5.70 on a lattice 16^ x 24, for each value of k, fitted pseudoscalar meson 
masses, time range of fit, and fit's x^ per degree of freedom obtained from propagators not using 
random sources. 



mass 



t range 



x'/d.o.f. 



0.1539 
0.1546 
0.1554 
0.1562 
0.1567 



0.4835(5) 
0.4456(5) 
0.3996(6) 
0.3496(7) 
0.3154(7) 



6- 9 
6- 9 
6- 9 
6- 9 

6- 8 



1.01 
0.86 
0.63 
0.38 
0.22 



TABLE VII. For (3 of 5.93 on a lattice 16^ x 14 x 20, for each value of k, fitted pseudoscalar 
meson masses, time range of fit, and fit's x^ per degree of freedom. 



mass 



t range 



x'/d.o.f. 



0.1539 
0.1554 
0.1567 



0.4820(4) 
0.3982(4) 
0.3147(4) 



8- 10 
8- 11 
8- 10 



0.40 
0.26 
0.10 



TABLE VIII. For (3 of 5.93 on a lattice 24*^, for each value of k, fitted pseudoscalar meson 
masses, time range of fit, and fit's x^ per degree of freedom. 
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mass 



t range 



xVd.o.f. 



0.1508 
0.1512 
0.1516 
0.1520 
0.1524 



0.3348(6) 
0.3094(6) 
0.2826(6) 
0.2541(6) 
0.2229(7) 



5- 14 
5- 14 
5- 14 
5- 14 
5- 14 



1.05 
1.22 
1.39 
1.56 
1.69 



TABLE IX. For /3 of 6.17 on a lattice 24^ x 20 x 32, for each value of k, fitted pseudoscalar 
meson masses, time range of fit, and fit's x^ per degree of freedom. 



mass 



t range 



x'/d.o.f. 



0.1485 
0.1488 
0.1491 
0.1494 
0.1497 



0.2564(6) 
0.2354(6) 
0.2133(7) 
0.1893(7) 
0.1630(8) 



13- 16 

13- 16 

14- 16 
14- 17 
14- 17 



0.96 
0.88 
0.00 
0.04 
0.12 



TABLE X. For /3 of 6.40 on a lattice 32^ x 28 x 40, for each value of k, fitted pseudoscalar 
meson masses, time range of fit, and fit's x^ P^r degree of freedom. 



mass 



t range 



x'/d.o.f. 



0.1600 
0.1613 
0.1625 
0.1638 
0.1650 



1.343(14) 
1.316(14) 
1.298(18) 
1.295(13) 
1.293(12) 



3- 6 
3- 5 
3- 5 
2- 4 

2- 4 



0.18 
0.55 
1.16 
0.00 
3.16 



TABLE XL For (3 of 5.70 on a lattice 12^ x 10 x 24, for each value of k, fitted scalar meson 
masses, time range of fit, and fit's x^ per degree of freedom. 



mass 



t range 



x'/d.o.f. 



0.1625 
0.1650 



1.299(11) 

1.287(12) 



3 

2 



5 

4 



0.32 
0.00 



TABLE XIL For (3 of 5.70 on a lattice 16^ x 24, for each value of k, fitted scalar meson masses, 
time range of fit, and fit's x^ psr degree of freedom. 
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mass 



t range 



xVd.o.f. 



0.1600 
0.1650 



1.325(15) 

1.278(3) 



5- 9 
2- 4 



0.19 
0.08 



TABLE XIII. For (3 of 5.70 on a lattice 16^ x 24, for each value of k, fitted scalar meson 
masses, time range of fit, and fit's x^ P^r degree of freedom obtained from propagators not using 
random sources. 



mass 



t range 



x'/d.o.f. 



0.1539 
0.1546 
0.1554 
0.1562 
0.1567 



0.860(4) 
0.837(4) 
0.815(5) 
0.812(6) 
0.818(6) 



4 
4 
4 
3 
3 



8 
8 
8 
7 
5 



2.39 
1.93 

1.48 
1.05 
0.84 



TABLE XIV. For /3 of 5.93 on a lattice 16^ x 14 x 20, for each value of k, fitted scalar meson 
masses, time range of fit, and fit's x^ per degree of freedom. 



K 



mass 



t range 



x'/d.o.f. 



0.1539 

0.1554 
0.1567 



0.851(9) 
0.806(4) 
0.779(6) 



7- 
4- 
4 



11 
11 

7 



0.33 

1.40 
1.48 



TABLE XV. For (3 of 5.93 on a lattice 24^, for each value of k, fitted scalar meson masses, 
time range of fit, and fit's x^ P^r degree of freedom. 



mass 



t range 



x'/d.o.f. 



0.1508 
0.1512 
0.1516 
0.1520 
0.1524 



0.574(4) 
0.559(5) 
0.546(6) 
0.538(8) 
0.547(7) 



6- 9 
6- 9 
6- 11 
6- 10 

4- 8 



0.25 
0.22 
0.22 
0.08 
0.26 



TABLE XVI. For (3 of 6.17 on a lattice 24^ x 20 x 32, for each value of k, fitted scalar meson 
masses, time range of fit, and fit's x^ per degree of freedom. 
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mass 



t range 



xVd.o.f. 



0.1485 
0.1488 
0.1491 
0.1494 
0.1497 



0.424(3) 
0.411(3) 
0.404(3) 
0.397(4) 
0.407(5) 



7- 13 
7- 18 
7- 18 
6- 13 
5- 8 



0.91 
0.86 
1.23 
1.03 
1.34 



TABLE XVII. For /? of 6.40 on a lattice 32^ x 28 x 40, for each value of k, fitted scalar meson 
masses, time range of fit, and fit's x^ per degree of freedom. 



(3 



lattice 



5.70 
5.70 
5.93 
5.93 

6.17 
6.40 



/3 



12^ X 10 X 24 
122 X 10 X 24 
16^ X 14 X 20 

24^ 
24^ X 20 X 32 
32^ X 28 X 40 



0.164382(23) 
0.164392(6) 
0.156391(11) 
0.156384(7) 
0.152167(11) 
0.149490(6) 



0.169538(70) 
0.169652(86) 
0.159062(15) 
0.159079(17) 
0.153833(18) 
0.150628(17) 



TABLE XVIII. Hopping constant at the strange quark mass and at zero quark mass. 



lattice 



mass 



5.70 
5.70 
5.93 
5.93 
6.17 
6.40 



(3 



12^ X 10 X 24 

16^ X 24 
16^ X 14 X 20 

24^ 
242 X 20 X 32 
322 X 28 X 40 



1.298(14) 
1.283(5) 
0.811(6) 
0.784(5) 
0.545(6) 
0.400(3) 



TABLE XIX. Scalar quarkonium mass interpolated to the strange quark mass. 



lattice 



mass 



t range 



xVd.o.f. 



lattice 



mass 



5.70 


122 X 10 X 24 


0.945(91) 


2- 


-4 


0.05 


16^ X 24 


0.955(15) 


5.93 


162 X 14 X 20 


0.788(21) 


1 - 


-4 


0.53 


16^ X 24 


0.781(11) 


5.93 


24^ 


0.774(23) 


1 - 


-4 


1.06 


16^ X 24 


0.781(11) 


6.17 


242 X 20 X 32 


0.577(39) 


2- 


-4 


0.00 


242 X 20 X 32 


0.559(17) 


6.40 


322 X 28 X 40 


0.397(35) 


4- 


- 6 


0.88 


322 X 30 X 40 


0.432(8) 



TABLE XX. For each /3 and lattice, scalar glueball mass, time range of fit, and fit's x^ P^i" 
degree of freedom, compared with masses obtained elsewhere from larger ensembles for the same 
P and nearly equal lattice sizes. 
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K mixing energy t range ^ jd.o.f . 

0.1600 0.167(15) ^T\ 0.43 

0.1613 0.180(14) 1 - 4 0.39 

0.1625 0.193(15) 1 - 4 0.37 

0.1638 0.205(14) 1 - 4 0.38 

TABLE XXI. For /3 of 5.70 on a lattice 12^ x 10 x 24, for each value of «;, fitted mixing energy, 

time range of fit, and fit's x^ per degree of freedom. 



K mixing energy t range X /'^•c/- 

0.1539 0.083(10) 2^5 0.33 

0.1546 0.088(10) 2 - 5 0.35 

0.1554 0.094(10) 2 - 5 0.40 

0.1562 0.099(10) 2 - 5 0.48 

0.1567 0.104(11) 2-5 0.52 

TABLE XXIL For /? of 5.93 on a lattice 16^ x 14 x 20, for each value of k, fitted mixing energy, 

time range of fit, and fit's x^ per degree of freedom. 



K mixing energy t range x /d-o.f. 

0.1539 0.105(19) ^n 0.77 

0.1554 0.115(17) 2-4 0.69 

0.1567 0.126(18) 2 - 4 0.52 

TABLE XXIII. For (5 of 5.93 on a lattice 24^, for each value of k, fitted mixing energy, time 
range of fit, and fit's x^ per degree of freedom. 



K, mixing energy t range x l^-O-f- 

0.1508 0.048(9) 3^5 0.76 

0.1512 0.051(9) 3 - 5 0.80 

0.1516 0.054(8) 3 - 5 0.85 

0.1520 0.057(8) 3 - 5 0.93 

0.1524 0.059(9) 3 - 5 1.08 

TABLE XXIV. For (3 of 6.17 on a lattice 24^ x 20 x 32, for each value of k, fitted mixing 
energy, time range of fit, and fit's x^ per degree of freedom. 
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mixing energy 



t range 



xVd.o.f. 



0.1491 

0.1494 
0.1497 



0.033(4) 
0.036(4) 
0.039(5) 



2- 6 
2- 6 
2- 6 



0.61 
0.59 
0.63 



TABLE XXV. For /3 of 6.40 on a lattice 32^ x 28 x 40, for each value of k, fitted mixing energy, 
time range of fit, and fit's x^ psr degree of freedom. 



/3 


lattice 


Eif^s) 


Eifin) 


E{fin)/E{fis) 


5.70 


12^ X 10 X 24 


0.211(16) 


0.258(19) 


1.22(3) 


5.93 


16^ X 14 X 20 


0.101(11) 


0.120(11) 


1.18(3) 


5.93 


24^ 


0.123(18) 


0.142(21) 


1.15(5) 


6.17 


242 X 20 X 32 


0.058(9) 


0.069(10) 


1.20(8) 


6.40 


322 X 28 X 40 


0.037(4) 


0.046(5) 


1.25(6) 



TABLE XXVI. Quarkonium-glueball mixing energy in lattice units for each /3 and lattice. 
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FIG. 1. Pseudoscalar effective masses and fitted mass for /? of 5.93 and k of 0.1554 on a lattice 
16^ X 14 X 20. 
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FIG. 2. Pseudoscalar effective masses and fitted mass for (3 of 5.93 and k of 0.1562 on a lattice 
16^ X 14 X 20. 
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FIG. 3. Pseudoscalar effective masses and fitted mass for (3 of 5.93 and k of 0.1567 on a lattice 
16^ X 14 X 20. 
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FIG. 4. Pseudoscalar effective masses and fitted mass for (3 of 5.93 and k of 0.1554 on a lattice 
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FIG. 5. Pseudoscalar effective masses and fitted mass for (3 of 5.93 and k of 0.1567 on a lattice 
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FIG. 6. Pseudoscalar effective masses and fitted mass for (3 of 6.40 and k of 0.1491 on a lattice 
32^ X 28 X 40. 
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FIG. 7. Pseudoscalar effective masses and fitted mass for (3 of 6.40 and k of 0.1494 on a lattice 
32^ X 28 X 40. 
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FIG. 8. Pseudoscalar effective masses and fitted mass for (3 of 6.40 and k of 0.1497 on a lattice 
32^ X 28 X 40. 
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FIG. 9. Scalar effective masses and fitted mass for (3 of 5.93 and k, of 0.1554 on a lattice 
16^ X 14 X 20. 
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FIG. 10. Scalar effective masses and fitted mass for (3 of 5.93 and k of 0.1562 on a lattice 
16^ X 14 X 20. 
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FIG. 11. Scalar effective masses and fitted mass for (3 of 5.93 and k of 0.1567 on a lattice 
16^ X 14 X 20. 
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FIG. 12. Scalar effective masses and fitted mass for (3 of 5.93 and k of 0.1554 on a lattice 24^. 
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FIG. 13. Scalar effective masses and fitted mass for (3 of 5.93 and k of 0.1567 on a lattice 24^. 
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FIG. 14. Scalar effective masses and fitted mass for (3 of 6.40 and k of 0.1491 on a lattice 
32^ X 28 X 40. 
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FIG. 15. Scalar effective masses and fitted mass for (3 of 6.40 and k of 0.1494 on a lattice 
32^ X 28 X 40. 
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FIG. 16. Scalar effective masses and fitted mass for (3 of 6.40 and k of 0.1497 on a lattice 
32^ X 28 X 40. 
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FIG. 17. Pseudoscalar quarkonium mass squared as a function of 1/k for f3 of 5.70 on a lattice 
12^ X 10 X 24. Results for 16^ x 24 are nearly identical. 
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FIG. 18. Pseudoscalar quarkonium mass squared as a function of 1/k for f3 of 5.93 on a lattice 
16^ X 14 X 20. Results for 24"^ are nearly identical. 
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FIG. 19. Pseudoscalar quarkonium mass squared as a function of 1/k for /3 of 6.17 on the lattice 
242 X 20 X 32. 
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FIG. 20. Pseudoscalar quarkonium mass squared as a function of 1/k for f3 of 6.4 on the lattice 
32^ X 28 X 40. 
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FIG. 21. Scalar quarkonium mass as a function of quark mass for [3 of 5.93. 
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FIG. 22. Scalar quarkonium mass as a function of quark mass for f3 of 5.7. 
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FIG. 23. Scalar quarkonium mass as a function of quark mass for /? of 6.17. 
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FIG. 24. Scalar quarkonium mass as a function of quark mass for f3 of 6.4. 
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FIG. 25. Lattice spacing dependence and continuum limit of the scalar ss mass, continuum 
limit of the scalar glueball masses, and one sigma upper and lower bounds on observed masses. 
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FIG. 26. Scalar glueball effective masses and fitted mass for /3 of 5.70 on a lattice 12^ x 10 x 24. 
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FIG. 27. Scalar glueball effective masses and fitted mass for /3 of 5.93 on a lattice 16^ x 14 x 20. 
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FIG. 28. Scalar glueball effective masses and fitted mass for (3 of 5.93 on a lattice 24^. 
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FIG. 29. Scalar glueball effective masses and fitted mass for f3 of 6.17 on a lattice 24^ x 20 x 32. 
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FIG. 30. Scalar glueball effective masses and fitted mass for (3 of 6.40 on a lattice 32^ x 28 x 40. 
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FIG. 31. Effective mixing energy and fitted mixing energy for /? of 5.93 and k of 0.1554 on a 
lattice 16^ x 14 x 20. 
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FIG. 32. Effective mixing energy and fitted mixing energy for /? of 5.93 and k of 0.1562 on a 
lattice 16^ x 14 x 20. 
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FIG. 33. Effective mixing energy and fitted mixing energy for /? of 5.93 and k of 0.1567 on a 
lattice 16^ x 14 x 20. 
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FIG. 34. Effective mixing energy and fitted mixing energy for /? of 5.93 and k of 0.1554 on a 
lattice 24''. 
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FIG. 35. Effective mixing energy and fitted mixing energy for /? of 5.93 and k of 0.1567 on a 
lattice 24''. 
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FIG. 36. Effective mixing energy and fitted mixing energy for fi of 6.40 and k, of 0.1491 on a 
lattice 32^ x 28 x 40. 
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FIG. 37. Effective mixing energy and fitted mixing energy for /? of 6.40 and k of 0.1494 on a 
lattice 32^ x 28 x 40. 
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FIG. 38. Effective mixing energy and fitted mixing energy for /? of 6.40 and k of 0.1497 on a 
lattice 32^ x 28 x 40. 
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FIG. 39. Glueball-quarkonium mixing energy as a function of quark mass for [3 of 5.93. 
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FIG. 40. Glueball-quarkonium mixing energy as a function of quark mass for [3 of 5.70. 
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FIG. 41. Glueball-quarkonium mixing energy as a function of quark mass for [3 of 6.17. 
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FIG. 42. Glueball-quarkonium mixing energy as a function of quark mass for [3 of 6.40. 
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FIG. 43. Lattice spacing dependence and continuum limit of the glueball-quarkonium mixing 
energy Ei^^Xg) and of the ratio E[nn)/E{iJ.s). 
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